Uncorrelated Random Networks 
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We define a statistical ensemble of non-degenerate graphs, i.e. graphs without multiple- and self- 
connections between nodes. The node degree distribution is arbitrary, but the nodes are assumed 
to be uncorrelated. This completes our earlier publication 1], where trees and degenerate graphs 
were considered. An efficient algorithm generating non-degenerate graphs is constructed. The 
corresponding computer code is available on request. Finite-size effects in scale-free graphs, i.e. 
those where the tail of the degree distribution falls like , are carefully studied. We find that in the 
absence of dynamical inter-node correlations the degree distribution is cut at a degree value scaling 
like A^''', with 7 = min[l/2, f/(/3 — 1)], where N is the total number of nodes. The consequence 
is that, independently of any specific model, the inter-node correlations seem to be a necessary 
ingredient of the physics of scale-free networks observed in nature. 
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I. INTRODUCTION 

This paper is a direct continuation of ref. The 
importance of defining statistical ensembles of random 
graphs in order to understand the geometry of wide 
classes of networks independently of any specific model 
was emphasized there. Concepts borrowed from field 
theory were used to define the ensemble of uncorrelated 
graphs and an algorithm generating such graphs was pro- 
posed. The general philosophy of our approach was illus- 
trated by focusing on a graph sub-ensemble, namely on 
an ensemble of connected trees with a scale-free degree 
distribution, where a number of hopefully interesting an- 
alytic results could be presented. But it should have been 
obvious that the adopted framework is of much broader 
applicability. Actually, we have explicitly stated that our 
algorithm generates efficiently not only trees, but also 
so-called pseudographs, called degenerate graphs in 0. 
However, we have also indicated that we have encoun- 
tered problems dealing with simple (i.e. non-degenerate) 
graphs. Hence the discussion of these non-degenerate 
graphs was postponed. We are now returning to the 
problem of defining the statistical ensembles of networks, 
which in the meantime has attracted the attention of 
other researchers 0-0 (at this point is it fair to men- 
tion also ref. , an early paper on uncorrelated graphs) . 

Although some overlap with is unavoidable, if this 
paper is to be self-contained, we would like to reduce the 
overlap to a minimum. The reader is invited to consult 
when he finds the discussion of this paper too sketchy. 

We shall not dwell much in introducing the subject. 
Let us recall that a graph is just a collection of nodes 
(vertices) connected by links (edges). It is a mathemat- 
ical idealization representing the various networks one 
encounters in nature, in social life, in engineering, etc. 
Quite often the pattern of connections between nodes 
looks fairly random. The concept of a random graph 
emerges quite naturally. 

Random graphs are interesting in themselves. There 



exists a classical theory of random graphs, a beautiful 
piece of pure mathematics ,7]. It turns out, however, that 
most large graphs one encounters in applications are not 
covered by this theory. The access - relatively recent - to 
the corresponding data triggered a rather intense activity 
(see the refs. ;9j] and references therein). 

Networks are also interesting considered from a 
broader perspective: it is useful to represent the archi- 
tecture, so to say the skeleton of many complex systems 
by an appropriate network. Hence, graphs are in a sense 
a gateway to the theory of complex systems, an exciting 
and promising new direction of research. 

The plan of this paper is as follows: In Sect. II we 
recall the definition of the ensemble of uncorrelated ran- 
dom graphs and we discuss the algorithm generating the 
graphs. In Sect. Ill we present the results of a sample 
of computer simulations, aimed to help understanding 
finite-size effects. The latter play a very important role 
as soon as the degree distribution has a fat tail. We 
explain the behavior of the data using an analytic ar- 
gument. We conclude in Sect. IV. For definiteness, we 
consider undirected graphs only, as in ilj. 



II. DEFINING THE ENSEMBLE 

Let us recall the construction proposed in ref. H . The 
partition function for the ensemble of random graphs is 
written as the formal integral defining a minifield theory 
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It will be seen that the non-negative constants p„ cor- 
respond to the degree distribution, while the auxiliary 
constants g, A and k control the dependence of Z on the 
number of nodes, links and loops, repectively. The inte- 
gral does not exist, but Z can be treated formally as a 
generating function in the Gaussian perturbation theory. 
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The main idea is to expand the exponential under the 
integral in ^ in powers of g : 

fdcbexp- [-(/)V2A] [1 + - V Pn^ 

K K ^ — ' 

•' n 

to get a series in g with well defined coefficients, viz. 
integrals with the Gaussian measure of integer powers 
of 4>. Each such integral is equal to a sum of contribu- 
tions, which can be represented graphically by the so- 
called Feynman diagrams ^3 • We have explained in ref. 
Q how such a diagram emerges, using a particular ex- 
ample. We do not have enough space to develop the 
point in more details. For those readers who are not 
conversant with field theory techniques we list the rules 
for constructing and calculating the Feynman diagrams 
corresponding to the term of order 0{g^) in jSJ: 

Each diagram has N labeled nodes. One should draw 
all topologically distinct diagrams, distributing degrees 
among nodes in all possible manners and connecting 
nodes pairwise. Self- and multiple-connections between 
nodes are allowed. Notice the similarity with the MoUoy- 
Reed construction With each diagram is associated 
a number, called the Feynman amplitude, determined by 
the following rules: each node of degree n contributes 
a factor ^Pnnl and each link contributes a factor kX. 
There is a symmetry factor 1/2 associated with every 
line connecting a node to itself and a symmetry factor 
1/m! associated with every m-tuple connection between 
nodes There is also a factor Zo/N\, the factorial 

being a remnant of the expansion of an exponential and 
Zq being the value of the Gaussian integral. 

Finally, the series representation of the partition func- 
tion reads 

N 1 ^ 

L,N ■ D ^ ' j=l 

where one sums over labeled diagrams D having a fixed 
numer of nodes and links, respectively N and L. S{D) 
is the product of factors 2 and factorials associated with 
self- and multiple-connections and rij is the degree of 
the j^^ node. One can show that the analogous series 
for log Z receives contributions of connected diagrams 
only. In this case the expansion in powers of k is a loop 
expansion: the leading term corresponds to L = -I- 1 
and comes from tree diagrams, the next term comes from 
one-loop diagrams etc. 

Our idea is to identify the Feynman diagrams of the 
toy model defined by eq. ^ with the graphs of a statis- 
tical ensemble. Indeed, Feynman diagrams are identical 
to graphs familiar to network community people, except 
that there are definite rules to calculate the correspond- 
ing weights. The minifield formulation enables one to 
summarize compactly the content of the model and has 



also the advantage of being a good starting point for an- 
alytical calculations, like those of ref. 

In the following, we always work with graph ensem- 
bles where N and L are fixed. Hence, up to an irrel- 
evant factor, the weight w of a labeled graph that is 
non-degenerate, i.e. such that nodes are neither multiply 
connected nor connected to themselves, is just 

N 

W'^Ylpn.rijl (4) 

In the presence of degeneracies one has to multiply the 
r.h.s. (right hand side) of by the factor S*"^ appearing 
in @. 

Eq. Q gives the weight of a microstate. Notice the 
factorized form and therefore the absence of non-trivial, 
dynamical correlations. Notice also that with the choice 
Pn oc (const)"/n! all non-degenerate graphs with the 
same N and L are equiprobable, because 

Y,n,^2L (5) 
j 

Thus, with a Poissonian p„ one recovers the classical 
graph ensemble of Erdos and Renyi. The ensemble un- 
der discussion is the most conservative generalization of 
the classical ensemble to the case of an arbitrary degree 
distribution. 

At this point the statistical ensemble is basically de- 
fined. However, in this paper, we wish to focus on 
non-degenerate graphs, which are the primary objects 
in graph theory. They correspond to a sub-ensemble of 
Feynman diagrams. In the conventional applications of 
field theory no specific recipe is formulated to single these 
diagrams out. Such a recipe is, however, needed here. 
Otherwise our definition of the ensemble would be too 
vague to be useful in applications. 

Before going farther let us outline the strategy we shall 
follow: as stated above our goal is now to complete the 
definition of the ensemble by the construction of an algo- 
rithm generating non-degenerate graphs. But we do not 
achieve this goal directly. First we construct, following 
ref m , an algorithm generating graphs that are degener- 
ate. Then, we show that the ensemble of these degenerate 
graphs is isomorphic, as far as the degree distribution is 
concerned, to the known model of balls- in-boxes [ij. Us- 
ing this mapping of one model on another we conclude 
that asymptotically the degree distribution _P„ in the en- 
semble of degenerate graphs is just Pn ■ Pn ^ Pn for 
N —>■ oo. Since we suspect that in this limit the de- 
gree distributions are the same for degenerate and non- 
degenerate graphs, we impose the appropriate constraint 
on the algorithm and perform a sample of computer ex- 
periments, to be described in the next section. The re- 
sults might seem surprising at first sight but a clear pic- 
ture eventually emerges when we estimate analytically, 
in the ensemble of degenerate graphs, the likelihood that 
a node is neither self- nor multiply-connected. 
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In a growing network model the construction of graphs 
is recursive and mimics a real physical process. In a 
static model like ours one does not refer to any physi- 
cal process. An ensemble is defined and the relative fre- 
quency of occurence of distinct graphs is fixed: If graphs 
A and B have weights P{A) and P{B), respectively, then 
they should be generated with a relative frequency equal 
to P{A)/P{B). Naively, one could imagine generating 
graphs uniformly in the space of graphs, accepting graph 
A, say, with probability P{A). However, such a uniform 
sampling is in practice very difficult to insure. Further- 
more, in an ensemble of very many graphs the accep- 
tance rate of the naive algorithm would be very small, 
since the normalized weight of any given graph is roughly 
speaking of the order of the inverse of the number of 
graphs. A clever idea is to introduce an appropriate 
random walk (Markov process) in the space of graphs, 
. . . — > Ak —> Ak+i Ak+2 which performs an 

importance sampling. The random walk is driven by the 
Markovian transition matrix P{A B). One can eas- 
ily show that if the transition matrix fulfills the detailed 
balance condition : 



P{A)P{A B)^ P{B)P{B A) 



(6) 



the frequency of the configuration Ak in the Markov pro- 
cess is proportional to P{Ak), provided one has moved 
away from the initial configuration. There are many 
P{A — » B) fulfilling the detailed balance condition for 
a given probability measure {P{A),yA}. One is free to 
choose any one. The simplest and popular choice 



P{A -> B)^ min{l, P{B) /P{A)} 



(7) 



is usually refered to as the Metropolis algorithm [l3| . 
The general idea of the method is problem independent. 
However, the choice of the proposed new configuration 
-B, given the current one A, is made by taking into ac- 
count the particularities of the problem at hand. Usually 
one proposes to change only slightly a small number of 
parameters in the current configuration. This insures a 
reasonable acceptance rate and minimizes the risk of per- 
forming time consuming calculation for nothing. 

The transition A ^ B logically involves two steps: 
one proposes B among all candidates and one accepts 
the proposal with a certain probability. One can write 
P{A — > i?) as a product of the probability Pc of choos- 
ing a particular candidate and of the probability Pa of 
accepting it : P[A ^ B) = Pc{A -> B)Pa{A B). 

Our algorithm works as follows. In the current con- 
figuration a random oriented link ij, the candidate for 
rewiring, is chosen. This is done with the probability 
1/2L. Then we select a vertex fc, with the probability 
1/N . The proposed move consists of rewiring ij into ik. 
Thus, P{A ^ B) = 1/2LN Pa{A B), and similarly 
for A ^ B. Inserting this into the detailed balance con- 
dition and dropping the factors 1/2L7V, identical on both 
sides of the equation, wc obtain 



which has the Metropolis solution for Pa{A B). Now, 
we use the fact that, according to Q), P{A) is a product 
of the node weights Furthermore, we observe that 



the rewiring changes rik ^ rik + 1 and nj — > nj - 
leaving the degrees of other nodes intact, to get 



1 only. 



pfA • (r^fc + l)!pn,-n jrij - l)!p„,-l 
Pa(A^B) = mm{l, ■ j } 

= mm{l,{nk + l)R{nk + l)/njR{nj)} (9) 



where R{n) — Pn/Pn-i- When rij = 1, the attempt 
is rejected, so that nodes with zero degree are never cre- 
ated. Notice, that we directly sample links to be rewired. 
The graphs produced by this algorithm are in general de- 
generate and multiply connected. It turns out, that the 
detailed balance condition and the way of sampling links 
insure that the symmetry factors in the weights of degen- 
erate graphs come out correctly. 

The presence of the factor (n^ -f l)/nj on the r.h.s. 
of means that the rewired nodes are sampled inde- 
pendently of their degree [T5l |. Furthermore, the rewiring 
depends on the node degrees only and is insensitive to the 
rest of the underlying graph structure. Hence, as far as 
the distribution of node degrees is concerned, the model 
is isomorphic to the well known balls-in-boxes model |14| , 
defined by the partition function 
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(10) 



and describing AI balls distributed with probability p„ 
among N boxes (in our case M = 2L). The constraint 
represented by the Kronecker delta on the r.h.s. of (|10|l is 
satisfied " for free" when cxo by virtue of Khintchin's 
law of large numbers, provided (n) = J2n ''^Pn/ J2nPn = 
M/N. The finiteness of (n) is always tacitly assumed in 
this paper. Hence, when the last condition is met the 
occupation number distribution of a single box P„ —^ pn 
for N ^ oo. 

Consequently, in the statistical ensemble including 
degenerate graphs the degree distribution becomes Pn 
asymptotically when the number of links is set to 



L = ^N{n) 



(11) 



PiA)PaiA -^B) = P{B)Pa{B ^ A) 



(8) 



When this condition is not met, the asymptotic degree 
distribution differs from p„, which is, in a sense, renor- 
malized. In particular, when L is smaller, this distribu- 
tion is Pn times an exponentially falling factor. When L 
is larger the situation depends on the shape of p„. When 
the latter is scale-free, p„ oc for large n, the distri- 
bution is (X Pn, except that an extra singular node with 
degree of order 0{N) shows up. These phenomena were 
discussed at length in the context of the balls-in-boxes 
model and also in our preceding work Q]. 

So far, only an algorithm generating degenerate graphs 
has been constructed. It is trivial to convert it into an 
algorithm producing non-degenerate graphs. It suffices 
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FIG. 1: The approach towards the hmiting rectangular shape 
of the calculated connectivity distribution P„: N = 100 
(squares) and 1000 (circles). The figure illustrates the claim 
that finite size correction fade away rapidly when the input 
connectivity distribution has a finite support. 
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FIG. 2: The evolution with A'^ of the average maximum de- 
gree of a graph for P — 2.1 (circles) and 3.0 (squares) and 
— 2^'^ to 2^*. This log-log plot illustrates the observation 
that for scale-free graphs with the exponent /3 approaching 
2, the maximum node degree evolves with A'^ at a rate close 
to that observed for /3 = 3, in variance with I12II . An ana- 
lytic argument explaining the observed behavior is given in 
the text. 



for that to add before the Metropolis test a few lines of 
code checking that the nodes i and k are neither identical 
nor linked. However, this check introduces a bias and it is 
not obvious what will be the degree distribution at finite 
N. 

A priori, the Metropolis test should insure that the 
number of nodes of degree n is close to Npn, provided 
the last number is large enough. And for fixed n it can be 
made arbitrarily large with a proper choice of N. Hence, 
a possible deviation of the degree distribution from p„ 
should be a finite size effect disappearing in the limit 
N —> oo when the couplings p„ are defined on a finite 
support. However, one has a problem when p„ has a fat 
tail. 

Let the degree distribution fall like n~'^, /3 > 2. For 
finite N it cannot fall like that indefinitely, there is a 
natural cut-off scaling as 



oc 



(12) 



The argument is well known: the expected number of 
nodes with n > Umax is less than unity. The presence of 
this cut-off was used by Dorogovtsev et al to explain 
why the observed scale-free networks are always charac- 
terized by a relatively small f3. 

Hence, coming back to the algorithm, there is always a 
range of n where fluctuations in the number of nodes are 
very large. Increasing TV does not help. Now, if certain 
fluctuations are systematically favoured by the constraint 
excluding degeneracies, then the resulting degree distri- 
bution can strongly deviate from the input weights p„. 
We dedicate a separate section to the discussion of this 
problem. 



III. FINITE-SIZE GRAPHS: DEGREE 
DISTRIBUTION 

Let us first consider a case where the support of p„ is 
finite, in order to check that in this case the problem is 
indeed under control. We perform a numerical experi- 
ment, setting for definiteness pn = 1/10 for n < 10 and 
Pn — otherwise, while L — 2.757V as dictated by Ullfl. 
The result is shown in Fig. 1: as expected, the conver- 
gence of the degree distribution towards the input one is 
very fast. 

Our next experiment is with scale-free graphs. As in 



Uur next exper 
ref. Q we set [l3 



^" = ^^- ^ V(/3-2)r(n + 2/3-3) " " ^''^ 

but the generated graphs are now non-degenerate with 
loops, instead of trees (the graphs are also, in general, not 
connected). Since (n) — 2, we also set L — N. We have 
chosen this example to illustrate a behavior which, as we 
shall argue in a moment, is generic. Since the same choice 
was made in ref. U , the reader can directly compare the 
results obtained in the two ensembles studied. 

The maximum node degree rimax has been measured in 
long runs. The results, for (3 = 2.1 and 3.0 are shown in 
Fig. 2. The important point is that the rate of increase 
of '^max with N is almost the same for the two values of 
(3, contrary to what one can read from the r.h.s. of H12|l. 
The exponent is slightly below the value 1/2, expected 
asymptotically (see below). The autocorrelation time in- 
creases roughly at the same rate, from about 860 (1810) 
to 5460 (10300) sweeps for /3 = 2.1 ( 3.0); in a sweep one 
attempts to rewire all the L links of the graph. The frac- 
tion of nodes and links belonging to the giant component 
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10° 10' 10' n 10' 

FIG. 3: The degree distribution for /3 = 2.1 and = 2^* 
(dot-dashed), 2^^ (dashed), 2^® (solid). The dotted line is the 
shape at infinite A'^. The cut-off scales like for /3 = 3 (see the 
next figure). Since the distribution is normalized to unity, 
this results in a deviation of its shape from the asymptotic 
one. 
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FIG. 4: The degree distribution for /3 = 3.0 and = 2" 
(dot-dashed), 2^^ (dashed), 2^® (solid). The dotted line is the 
shape at infinite A^. 



decreases very slowly from about 0.57 and 0.78 (0.68 and 
0.81) at TV = Q 55 j^j^j Q ( q gg q 

N = 2'\ 

The shape and the evolution with N of the degree dis- 
tribution is shown in Figs. 3 and 4. It is manifest that the 
distribution does approach the expected limit. However, 
the approach is very slow and non-uniform, especially for 
/? = 2.1 . 

The results of the computer experiment can be under- 
stood using the following heuristic arguments: 

Let n 3> 1 be the degree of a node in the tail, so that 
there is no more than one node of that degree. We con- 



sider the ensemble of degenerate graphs and we estimate 
the fraction of graphs where the node in question has 
no self- or multiple-connections. The counting is easily 
performed by considering the symmetric adjacency ma- 
trix Cij — #links ij (every undirected link contributes 
to two elements of C). Obviously, dj is or 1 for non- 
degenerate graphs, but can take any integer value for 
degenerate ones. Let m be the label of our hub. We 
count the adjacency matrices satisfying C'jm = n. 

We take the hmit N,L ^ 00, L/N fixed and n/N ^ 
(notice, that it is not assumed that n is kept fixed when 
N increases, it is only less than the cut-off given by eq. 
(|12|l ). In this limit, the number of degenerate graphs 
we are interested in is proportional to the number of 
ways to place n elements in N cells, possibly putting 
several elements in the same cell. The standard fomula 
of combinatorics does not apply because one has to take 
care of the contribution of the symmetry factors to graph 
weights. The number of graphs corrected by the appro- 
priate weights is 

1 1 ^ 

{fe.,} j i=l 

= _L / rf„e-ma+(iV-l/2)e- (^5) 

27r J 

where kj = Cjm ■ There is a self-connection when k.„i > 
and a multiple connection when, for j ^ m, kj > 1. The 
symmetry factors are those mentioned at the place where 
we summarize the Feynman rules. We have replaced the 
i5-function by its Fourier transform and performed the in- 
dependent summations over kj . Saddle-point integration 
yields 

WdegH cx — (16) 

A similar counting of graphs where our node is neither 
self- nor multiply connected yields 

/TV - 1\ TV" ^2 
W^non-dcgH = ^ j °'^''''P^"2]V^ ^^^^ 

The ratio is [3 

2 

Tl 

M^non-deg/W^deg « exp(- — ) (18) 

Hence, the entropy of the non-degenerate graphs is dra- 
matically reduced compared to that of graphs with de- 
generacies and this reduction depends strongly on n. 
This is the origin of the bias mentioned in the preced- 
ing section: 

Suppose we have generated a very large sample of 
scale-free graphs without any care for degeneracies. For 
each such graph the histogram of degrees contains a tail 
of sparsely located columns of unit height. However, 
we know, because of the mapping on the balls-in-boxes 
model, that the histogram for the whole sample has the 
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shape of Pn for n < Umax- Now, what happens when we 
exclude degeneracies, specifically when we check whether 
a node in the tail of the histogram does have the forbid- 
den connections? Eq. (|18|l tells us that the rejection rate 
is non-uniform in n and that nearly all candidate graphs 
will be rejected when n is large compared to N^^"^. 

Thus, to the exent one can neglect the fluctuations of 
graph weights, the degree distribution in non-degenerate 
graphs is expected to be cut by an additional factor, 
roughly behaving like exp (— n^/2iV). In the absence of 
dynamical correlations and for /3 < 3 the cut-off in non- 
degenerate graphs is expected to scale like iV^/^ and not 
like the r.h.s. of eq. (|12|l . i.e. independently of (3 ! Im- 
posing non-degeneracy one generates, at finite N , para- 
site correlations whose manifestation is the violation of 
eq. ()12|) . Notice, that this entropy argument does not ap- 
ply to trees, which within the full graph ensemble have 
a nearly vanishing entropy and which, as shown in Hll, 
have their own mapping on balls- in-boxes. 

The conclusion of this section is that our algorithm 
is efficient in generating all classes of graphs for any 
given degree distribution. In the specific, but most in- 
teresting case, of the maximally random non-degenerate 
scale- free graphs whose degree distribution falls like , 
with /3 < 3, the tail of the degree distribution is cut at 
f^max ~ N^/'^. This effect prevents the finite networks 
from fully developing the large part of the tail of their a 
priori expected degree distribution. This is a feature of 
the model, a result of the absence of dynamical correla- 
tions, and not a failure of the algorithm 19j. 

It is known that nodes are correlated in some growing 
network models with a scale- free behavior Here we 
find, independently of any specific model, that dynamical 
correlations seem to be needed for some networks with 
/3 « 2 to be observed in the real world. We also ob- 
serve that the thermodynamic limit N ^ oo oi the max- 
imum entropy graph model can be rather non-uniform 
and therefore somewhat tricky. 

IV. SUMMARY AND CONCLUSION 

There is much activity in producing new growing net- 
work models. Such models are invaluable for illustrating 
some basic dynamical mechanisms, like the preferential 
attachment rule |2ll |. However, in order to understand 
the generic geometries of wide classes of networks it is 
perhaps worthwhile to adopt a complementary approach, 
consistently defining the corresponding statistical ensem- 
bles and working with these static ensembles using the 
standard tools of equilibrium statistical mechanics. 

With this motivation in 1] we have defined a statis- 
tical ensemble for an arbitrary node degree distribution. 
This ensemble is the most random possible, but clearly 



not the most general: we assumed that node degrees are 
independent, to the extent that this is possible when the 
number of nodes and links is fixed. In ref. we have 
rapidly focused on trees, indicating, however, that the 
approach extends to more complicated graphs. The dis- 
cussion of the ensemble of non-degenerate graphs, those 
without multiple- and self-connections, has been post- 
poned. The present paper completes ref. 

In particular, we have constructed here an efficient al- 
gorithm, easily implemented on a computer, which en- 
ables one to generate non-degenerate random graphs. We 
are ready to share our computer code with interested peo- 
ple, on request. 

Another simple algorithm has been proposed in ref 
lllj l and used subsequently: one first generates from 
a given distribution p„ a set of node degrees {rij} and 
uses these numbers to construct auxiliary graphs, each 
with one node and rij half-links. In the second step one 
connects the half-links at random. The resulting graph is 
usually degenerate Imposing the absence of degen- 
eracies is in this approach very tedious, especially that a 
non-degenerate graph may just not exist for a given set 

We have studied in detail the behavior of the de- 
gree distribution of finite-size graphs. In the absence of 
dynamical inter- node correlations, for generic scale- free 
non-degenerate graphs (but not trees) this distribution is 
cut at 

rimax cx , 7 = min[l/2, l/(/3 - 1)] (19) 

at asymptotically large N . It appears that a fat tail with 
13 rather close to 2 observed in some data could hardly 
show up if dynamical inter-node correlations were absent. 
And indeed, non-trivial correlations are present in models 
of growing networks using the preferential attachment 
recipe 17] . 

It is certainly very important to develop a theory of 
correlated networks. Other authors 0, Q have very re- 
cently made interesting explicit proposals in that direc- 
tion. Our approach can also be rather easily generalized 
to include dynamical correlations. This does not mean 
that the physics of correlated networks is transparent to 
us at the present time. We shall hopefully return to the 
problem of correlations in another publication. 
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